one  FILE  COPY  M  A 1 253  36 


TECHNICAL  REPORT  ARBRL-TR-02465 


TWO-PHASE  VISCOUS  FLOW  MODELING  OF 
INTERIOR  BALLISTICS,  ALGORITHM,  AND 
NUMERICAL  PREDICTIONS  FOR  AN 
IDEALIZED  LAGRANGE  GUN 


James  A.  Schmitt 
Norman  E.  Banks 
Csaba  K.  Zoitani 
Thomas  L  Mann 
Howard  J.  Gibeling 

January  1983 


US  ARMY  ARMAMENT  RESEARCH  AM  DEVELOPMENT  COMMAND 
BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND,  MARYLAND 


Appro  v*d  public  itlNUi  distribution  unllattod. 


^3  02  014  052 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Secondary  distribution  of  this  report  is  prohibited. 


Additional  copies  of  this  report  may  be  obtained 
from  the  National  Technical  Information  Service, 

U.  S.  Department  of  Commerce,  Springfield,  Virginia 

22161. 


The  findings  in  this  report  are  not  to  be  construed  as 
an  official  Department  of  the  Array  position,  unless 
so  designated  by  other  authorized  documents. 


The  u*a  of  trad*  norm*  or  mjnufaoturar* '  nan**  in  tki*  report- 
do**  not  constitute  inJoremnent  of  any  oommraial  product . 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  or  THIS  PAGE  (When  Data  Bnlr rod) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

i.  rIpoKt  mumper  a.  govt  accession  no. 

Technical  Report  ARBRL-  TR-02465  ,$£)-/)  /if 

».  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (md  Submit) 

Two- Phase  Viscous  Flow  Modeling  of  Interior 
Ballistics,  Algorithm,  and  Numerical  predictions 
for  an  Idealized  Lagrange  Gun 

8.  TYPE  OP  REPORT  k  PERIOD  COVERED 

Technical  Report 

8-  PERFORMING  ORG.  REPORT  NUMBER 

7.  author^ 

James  A.  Schmitt,  Norman  E,  Banks,  Csaba  K. 

Zoltani,  Howard  J.  Gibeling*,  Thomas  L,  Mann 

B.  CONTRACT  OR  GRANT  NUMBER^.) 

•  ■  PERFORMING  ORGANIZATION  NAME  AND  ADORES-. 

US  Army  Ballistic  Research  Laboratory 

ATTN:  DRDAR-BLI 

Aberdeen  Proving  Ground,  MD  21005 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  8  WORK  UNIT  NUMBERS 

1L161 102AH43 

It.  CONTROLLING  OrriCE  NAME  ANO  ADDRESS 

US.  Army  Armament  Research  and  Development.  Command 

US  Army  Ballistic  Research  Laboratory  (DRDAR-BL) 
Aberdeen  Proving  Ground,  MD  21005 

11.  REPORT  DATE 

January  1983 

IS.  NUMBER  OP  PAOES 

38 

U.  MONITORING  AGENCY  NAME  •  ADDRESS?/?  dlltotont  from  Cantroltlnl  Otlloe) 

IS.  SECURITY  CLASS,  (ol  I bit  report) 

UNCLASSIFIED 

(s«.  OECL  ASSI  PIC  ATION/  DOWN  GRADING 
SCHEDULE 

u.  distribution  statement  (oi  mi*  Report; 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (ol  tfio  atl.tr act  antarad  In  Block  20,  II  d» /(trail  from  Report) 

IS.  SUPPLEMENTARY  NOTES 

•Scientific  Research  Associates,  Inc. 

Glastonbury,  CT  06033 

IS.  KEY  WORDS  (Continue  on  nntM  eldo  It  neeeeeory  and  identify  by  block  number) 

Interior  Ballistics  ALPHA  Adiabatic  Wall 

Lagrange  Gun  Laminar 

Two-phase  Turbulence 

Viscous  Flow  Isothermal  Wall 

SO.  ABSTRACT  (CmtBw  om  eoeoroo  • Mr  W  imma,  and  idomtlty  by  block  number)  .»  j  j 

A  new  state  of  the  art  algorithm,  ALPHA,  for  the  simulation  of  the 
multiphase,  multidimensional,  unsteady,  compressible,  viscous,  non-reactive, 
interior  ballistics  flow  in  a  gun  tube  behind  an  accelerating  projectile  is 
described.  The  paper  contains  discussions  of  the  physical  processes  in  a 
real  gun  environment,  of  the  mathematical  model  of  these  phenomena,  and  of 
the  numerical  technique  for  solving  the  equations.  The  algorithm  allows  the 
inclusion  of  several  submodels,  such  as,  heat  transfer  and  (CONTINUED) 

MMrt.EfiHTgTgn _ 

SECURITY  CLASSIFICATION  or  THIS  RASE  f Btwm  Data  Bnlorod) 


DD  |  1473  cornoM  or  t  mov  m  is 


■  TE 


_ UNCLASSIFIED _ 

SZCURITY  CLASSIFICATION  Of  THU  FAOKWRl  Dm*  MmHroHj, 

20.  ABSTRACT  (Continued) « 

turbulence,  This  permits  the  determination  of  the  effects  of  these 
submodels  on  the  flow.  Numerical  results  of  an  idealized,  one-chase 
ballistic  cycle  are  given.  Some  of  the  significant  results  include  the 
existence  of  a  concentrated  region  of  high  temperature  near  the  juncture 
of  the  projectile  base  and  tube  wall,  the  increase  of  the  displacement 
thickness  by  a  factor  of  at  least  three  over  most  of  the  tube's  length 
when  turbulence  effects  are  included,  and  the  degradation  of  the 
projectile  velocity  by  approximately  ten  percent  under  an  isothermal 
cold  wall  condition. 


i 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  FAO Date  MnltnH) 


TABLE  OF  CONTENTS 


Page 

LIST  OF  FIGURES .  5 

I.  INTRODUCTION .  7 

II.  GOVERNING  EQUATIONS .  9 

III.  ALGORITHM  DESCRIPTION . 12 

IV.  THE  MODEL  PROBLEM . 15 

A.  Laminar  Flow-Adiabatic  Wall  Simulation . 16 

B.  Turbulent  Flow-Adiabatic  Wall  Simulation . 20 

C.  Laminar  Flow-Isothermal  (Cold  Wall)  Simulation . 22 

D.  Comparison  Among  the  Simulations. . . . ...23 

V.  CONCLUSIONS  AND  FUTURE  EFFORTS . 24 

REFERENCES . 26 

APPENDIX . 29 

NOMENCLATURE . 35 

DISTRIBUTION  LIST . 37 


t  • ■  r-  n  nt  on  Fc 

"  G!U*I 
'  TAB 
■  ».  runic  "d 

. :  f  leatio 


i  !■ .  ribut Ion/ 

I  •  vatlabi lity  Codss 
j  Avail  and/or 

Spec  1 «1 


3 


COPY 

INSPECTED 


LIST  OF  FIGURES 


Illustration  of  the  Ballistic  Cycle. 

( Adaptation  from  Ref .  1)........... . . 

Pressure  Histories  at  the  Breech  (  — )  and  at  Projectile 

Base  ( - )  and  Their  Comparisons  with  Analytically 

Predicted  Initial  Arrival  Times  of  the  Rarefaction 
Wave . . . 

Velocity  Boundary  Layer  Profiles  for  the  Laminar 

Flow-Adiabatic  Wall  Calculation . 

Thermal  Boundary  Layer  Profiles  for  the  Laminar 

Flow-Adiabatic  Wall  Calculation . 

Velocity  Boundary  Layer  Profiles  for  the  Turbulent 

Flow-Adiabatic  Wall  Calculation . 

Thermal  Boundary  Layer  Profiles  for  the  Turbulent 

Flow-Adiabatic  Wall  Calculation . . . 

Velocity  Boundary  Layer  Profiles  for  the  Laminar 

Flow-Isothermal  (Cold  Wall)  Calculation . 

Thermal  Boundary  Layer  Profiles  for  the  Laminar 

Flow- isothermal  (Cold  Wall)  Calculation . 

Comparisons  of  the  Displacement  tfiicknesses  for  the  Three 
Simulations  of  the  Model  Lagrange  Gun  Problem  at 
t*2.4  ms . . . 

Differences  in  Projectile  Velocities  as  a  Function  of 
Time,  Flow,  and  Boundary  Conditions . 


PREVIOUS  PAGE 
IS  BLANK 


I .  INTRODUCTION 


To  advance  the  state-of-the-art  of  gun  system  development,  especially 
those  with  high  muzzle  velocity,  a  detailed  understanding  of  the  phenomena 
occurring  within  the  gun  tube  is  needed.  Current  theoretical  and  experimental 
research  at  the  Interior  Ballistics  Division  of  the  Ballistic  Research 
Laboratory  and  Scientific  Research  Associates,  Inc.,  addresses  this  problem 
with  special  emphasis  on  the  modeling  of  the  unsteady,  multiphase  aspect  of 
the  interior  ballistic  flows.  The  ballistic  cycle,  as  illustrated  in 
Figure  1 , 1  commences  with  the  ignition  of  the  propellant  charge  and  terminates 
with  the  exit  of  the  projectile  and  the  emptying  of  the  hot  gases  and  any 
unburned  propellant  from  the  gun  tube.  Typically,  an  ignition  of  the 
propellant  bed  is  started  by  a  hot,  gas-particle  flow  from  the  igniter.  As 
the  flame  spreads  throughout  the  propellant  bed,  gases  are  generated,  pressure 
waves  evolve,  fluidization  of  the  propellant  bed  begins,  and  subsequently, 
this  multiphase  flow  proceeds  down  the  tube  behind  the  nonuniformly 
accelerated  projectile.  Once  there  is  appreciable  axial  projectile 
displacement,  several  ancillary  phenomena  can  be  observed.  These  include 
projectile/tube  interaction  manifested  by  balloting  and  tube  vibration  as  well 
as  the  leakage  of  propellant  gases  around  and  ahead  of  the  moving  projectile. 
These  processes  are  not  considered  in  this  report. 
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Figure  1.  Illustration  of  the  Ballistic  Cycle 
(Adaptation  from  Ref.  1) 
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Our  interest  is  focused  on  the  flow  within  a  gun  bounded  by  the  breech, 
the  gun  tube  wall,  and  the  projectile  base.  From  a  global  point  of  view, 
interior  ballistics  phenomena  can  be  divided  into  two  coupled  parts:  the  two- 
phase  flow  behind  the  projectile  and  the  motion  of  the  projectile.  The  two- 
phase  flow  provides  a  force  which  partially  determines  the  motion  of  the 
projectile,  and  the  projectile  base  is  an  accelerating  boundary  for  the  two- 
phase  flow  field.  From  our  standpoint,  the  projectile  motion  can  be  treated 
by  a  lumped  parameter  analysis  in  which  the  projectile  moves  axially  with  a 
possible  fixed  rotation,  if  in  a  rifled  gun  tube,  and  with  a  perfect  seal  and 
no  gas  leakage.  Applying  Newton's  second  law  to  the  projectile  axial  motion, 
we  obtain  an  ordinary  differential  equation.  The  forces  on  the  projectile 
include  the  area  integral  of  the  pressure  at  the  base  of  the  projectile  plus 
the  wall  friction  force,  the  force  needed  to  engrave  the  rotating  bands  on  the 
grooves  of  a  rifled  gun  tube,  the  rotational  force,  and  the  air  resistance. 

The  last  four  forces  are  usually  prescribed  for  a  given  gun  and  projectile. 

Until  recently,  the  most  sophisticated  modeling  of  major  portions  of  the 
ballistic  cycle  has  been  limited  to  quasi-one-dimensional,  inviscid,  two-phase 
flow  analyses2-*  .  By  assuming  a  cylindrically  symmetric  flow,  parts  of  the 
ballistic  cycle  can  now  be  modeled  two-dimensionally .  Gough8  has  recently 
developed  an  inviscid,  two-dimensional  model  to  study  these  phenomena. 

ALPHA7'8  is  the  first  two-phase,  two-dimensional,  Navier-Stokes  model 
developed  to  simulate  the  ballistic  cycle.  None  of  the  above  models  includes 
chemical  effects. 


2P.S.  Gough  and  F.J.  Zwarte,  "Modeling  Heterogeneous  Two-Phase  Reacting 
Flow."  AIM  Journal.  Vol .  17,  Ho.  1,  pp.  17-25,  1979. 

ZK.K.  Kuo,  J.H.  Koo ,  T.R.  Davis,  and  G.R.  Coatee ,  " Transient  Combustion  in 
Mobile  Gas-Permeable  Propellants,"  Acta  Astronautica ,  Vol.  3,  pp.  573-591 , 
1976. 


4E.B.  Fisher,  K.W.  Graves,  and  A.P.  Trippe,  "Application  of  a  Flame  Spread 
Model  to  Design  Problems  in  the  155-mi  Propelling  Charge, "  1 2th  « TAHNAF 
Combustion  Meeting.  CPI A  Pub.  273,  Vol .  I,  pp.  199-219 ,  December  1975. 

5H.  Krier  and  S.S.  Gokhale,  "Modeling  of  Convective  Mode  Combustion  Through 
Granulated  Propellant  to  Predict  Detonation  Transition, "  AT A A  Journal, 

Vol.  16,  Ho.  2,  pp.  177-183,  1978. 

8P.S.  Gough,  "A  Two-Dimensional  Model  of  the  Interior  Ballistics  of  Bagged 
Artillery  Charges,"  USA  ARRADCOM/Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground,  MD,  ARBRL-CR-00452 ,  April  1981. 

7H.J.  Gibeling,  R.C.  Buggeln,  and  H.  McDonald,  "Development  of  a  Two- 
Dimensional  Implicit  Interior  Ballistics  Code,"  USA  ARRADCOM/Ballistic 
Research  Laboratory ,  Aberdeen  Proving  Ground ,  MD,  ARBRL-CR-00411,  January 
1980. 


O  # 

H.J.  Gibeling  and  H.  McDonald,  "Development  of  a  Two-Dimeneinal  Implicit 
Interior  Ballistics  Code,"  USA  ARRADCOM/Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  MD,  ARBRL-CR-00451,  March  1981. 
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The  purpose  of  this  report  is  to  present  a  two-phase,  viscous  model  and 
to  detail  the  calculations  for  an  idealized  one-phase  ballistic  cycle. 

Although  these  first  applications  of  ALPHA  exclude  the  second  phase,  important 
effects  of  the  basic  mechanisms  of  heat  transfer  to  the  gun  tube  wall  and  of  a 
turbulent  flow  are  isolated  and  resolved,  in  Section  II,  the  governing 
equations  which  provide  a  model  of  the  two-phase  interior  ballistics  flow  are 
listed.  The  description  of  the  numerical  scheme  to  solve  these  equations  is 
also  given.  The  results  of  three  simulations:  the  laminar  flow-adiabatic 
wall,  the  turbulent  flow-adiabatic  wall,  and  the  laminar  flow-isothermal  wall, 
are  then  presented.  Section  V  assays  the  results  and  discusses  future  work. 


II.  GOVERNING  EQUATIONS 

Because  of  the  complex  nature  of  the  two-phase  flow  behind  the 
projectile,  simplifying  assumptions  must  be  made  prudently  in  order  that 
significant  phenomena  are  not  neglected  in  the  model.  Therefore,  we  start 
with  a  general  formulation  and  introduce  simplifications  only  when 
warranted.  Due  to  the  large  number  of  propellant  grains,  of  the  order  of  a 
few  thousands,  the  mathematical  description  is  restricted  to  the  bulk  or 
average  properties  of  the  flow.  An  averaged  variable  represents  the  integral 
of  the  product  of  a  microscopic  variable  and  a  weighting  function,  which 
reflects  the  influence  of  remote  points  on  the  average  value,  over  time  and 
space.  Likewise,  the  equations  of  continuum  mechanics  are  replaced  with  a  set 
of  averaged  equations  in  the  averaged  variables  at  each  spatial  position  and 
time.  The  derivation  of  these  averaged  equations  is  discussed  in  Refs.  7  and 
8.  The  constitutive  laws  represent  relations  between  the  averaged  variables 
and  the  different  phases,  and  are  given  in  the  Appendix.  The  averaged, 
Navier-Stokes  equations  are  summarized  here  in  dimensional  form. 

The  gas  and  solid  phase  continuity  equations  are 
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respectively,  where  the  porosity  g  is  the  ratio  of  volume  occupied  by  the  gas 
phase  to  the  total  volume.  The  averaged  densities  and  velocities  of  the  gas 
and  solid  phases  are  denoted  by  p,  U,  p  ,  U  ,  respectively.  The  interphase 
mass  transfer  for  t*  e  gas,  T,  is  due  to^  th£  burning  propellant  and  is  given 
by  Eq.  (19).  (Eqs.  (19)  through  (40)  are  in  the  Appendix.)  The  gas  and  solid 
phase  momentum  equations  are 
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where  P  is  the  averaged  pressure.  The  tensors,  tt  and  tt  ,  Eqs.  (23)  and  (26), 
are  the  laminar  and  turbulent  stress  tensors,  respectively.  The  constitutive 
relations  fgr  grain  surface  area,  Sp,  and  volume,  V  ,  the  interphase  drag 
relation,  <F>,  and  the  intergranular  stress  relation,  R  ,  are  giv^n  by  59s. 
(20),  (21),  (29),  and  (30),  respectively.  The  velocity  vectors,  U  and  U  , 
have  three  components;  a  radial,  an  axisymmetric  angular,  and  an  axial.  ^ 
Axisymmetric  swirl  may  be  important  for  the  case  of  the  projectile  rotation. 
The  gas  phase  energy  equation  is 
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where  h  is  the  enthalpy  of  the  gas.  The  dissipation  function,  *,  the 
turbulent  kinetic  energy  dissipation  rate,  e,  the  energy  transfer  term  between 
the  solid  and  gas  phases,  A,  and  'hi  laminar  and  turbulent  heat  flux 
vectors  q  and  q  are  given  by  Eqs.  (31),  (28),  (32),  (34),  and  (35), 
respectively.  The  symbol  D/Dt  denotes  the  material  derivative.  The  Noble- 
Abel  equation  of  state,  (38),  is  used  for  the  gas.  The  turbulent  kinetic 
energy  is  given  by 
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where  a.  is  set  equal  to  one.  The  rate  of  strain  tensor  D  and  the  turbulent 
k  t 

viscosity  m  are  given  by  Eqs.  (24)  and  (27),  respectively. 

By  excluding  the  chemical  reactions  in  the  present  two-phase  flow 
analysis,  the  gas  phase  species  and  gasified  propellant  species  mass  fractions 
are  not  required.  However,  in  order  to  consider  the  effects  of  several  types 
of  propellents  within  the  tube,  transport  equations  for  the  reciprocal  of  the 
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gas  mixture  molecular  weight  M  and  specific  heat  at  constant  pressure  cp  are 
solved.  The  governing  equations  for  M  and  Cp  are 
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where  M  and  (Cp)p  are  the  reciprocal  molecular  weight  and  specific  heat  at 

constant  pressure  of  the  propellant,  respectively,  T  ■  y  __/Sc  __  , 

T  m  eff  erf 

Meff  “  y  +  M  '  and  Sceff  *  °*9, 


Presently,  we  assume  that  the  propellant  grains  are  spherical.  To 
determine  ignition,  the  surface  temperature  of  these  spheres  must  be 
calculated.  If  we  assume  that  the  penetration  depth  of  the  thermal  wave  into 
the  grain  is  small  compared  to  the  grain  diameter,  a  one-dimensional 
approximation  can  be  used  to  obtain  the  propellant  surface  temperature.  The 
appropriate  equation  for  the  propellant  temperature  Tp  is 
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where  r  is  the  radial  location  within  the  spherical  grain  and  dp  is  its 
thermal  diffusivity.  Once  the  surface  temperature  of  the  grain  exceeds  the 
ignition  temperature  of  the  propellant,  the  propellant  grain  is  assumed  to 
burn  until  it  is  completely  consumed.  The  actual  burning  is  modeled  as  a 
regression  of  the  surface  of  the  propellant  which  results  in  the  deposition  of 
mass,  momentum,  and  energy  into  the  gas.  The  equation  for  the  particle 
radius,  rp,  including  turbulent  diffusion,  is 
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where  the  regression  rate  <d>  is  given  by  Bq.  (22). 


The  effects  of  the  primer  which  deposits  hot  gases  into  the  propellant 
bed  to  begin  ignition  can  be  modeled  in  several  ways.  One  approach  is  to 
treat  the  "primer  source"  as  a  boundary  condition  or  initial  condition  either 
near  the  tube  centerline  for  a  center  core  igniter  or  near  the  breech  end  for 
a  base  igniter9. 
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The  two-phase  governing  equations,  Eqs.  (1)  through  (8)  and  (10),  are 
written  as  a  general  system  of  time-dependent  nonlinear  partial  differential 
equations.  The  d^perjdent  variables  are:  the  components  of  gas  and  solid 
phase  velocities  U,  Up,  partial  densities  ap,  (l-a)p ,  gas  phase  static 
enthalpy,  h,  specific  neat,  c  ,  reciprocal  mixture  molecular  weight,  M,  the 
turbulent  kinetic  energy,  k,  and  solid  phase  particle  radius,  r  .  Subsequent 
to  the  solution  of  this  coupled  system  of  gas-solid  phase  equations,  the  9olid 
particle  surface  temperature  is  determined  from  the  solution  of  Bq.  (9). 


III.  ALGORITHM  DESCRIPTION 

For  the  algorithm  an  implicit  formulation  is  favored  in  view  of  the 
desired  boundary  layer  resolution  along  the  gun  barrel.  To  treat  the 
nonlinearities  of  the  governing  equations,  '"aylor-series  expansions  in  time 
are  used  to  linearize  the  system  to  the  same  order  accuracy  as  the  temporal 
discretization.  To  efficiently  eliminate  the  large  block  banded  matrix,  a 
splitting  of  the  matrix  into  a  sequence  of  more  easily  treated  component 
matrices  is  performed.  Because  of  this  construction,  the  resulting  algorithm 
is  termed  a  split  _linearized  block  jjnplicit  (LBI)  scheme. 

The  numerical  algorithm  is  described  in  detail  in  Refs.  10  and  11. 
Following  Ref.  12,  a  brief  derivation  of  the  3cherr.e,  which  includes 
consideration  of  intermediate  boundary  conditions  required  by  split  schemes, 
is  now  given.  By  considering  these  intermediate  boundary  conditions,  a 
commutative  condition  on  the  boundary  conditions  is  identified  and  it  is  shown 
that  the  treatment  of  these  intermediate  boundary  conditions  can  affect 
transient  accuracy,  and  may  even  cause  errors  in  steady  solutions. 

To  Illustrate  the  derivation,  let 

($n+1  .  $°)/At  -  8D(Jn+1)  +  ( 1-0>D( )  +  0 [ At2 ,  ( 0-1 /2) At]  (11) 


approximate  3$/3t  *  D($),  a  system  of  £i me -dependent  nonlinear  partial 
differential  equations  for  the  vector  $  of  dependent  variables,  where  D  is  a 
multidimensional  vector  spatial  differential  operator,  and  t  is  a  time 


W.R.  Briley  and  H.  McDonald,  "Solution  of  the  Multidimensional  Compressible 
Navier-Stokes  Equations  by  a  Generalised  Implicit  Method ,  "  Journal  of 
Computational  Physics ,  Vol.  24,  No.  4,  pp.  372-397,  1977. 

^W.R.  Briley  and  H.  McDonald,  "On  the  Structure  and  Use  of  Linearised  Block 
ADI  and  Related  Schemes,"  Journal  of  Computational  Physics,  Vol.  34,  No.  1, 
pp.  5 4-73,  1980. 

5> 

JH.  McDonald  and  W.R.  Briley,  "Some  Observations  on  Numerical  Solutions  of 
the  Three-Dimensional  Navier-Stokes  Equations, "  presented  at  Symposiien  on 
Numerical  and  Physical  Aspects  of  Aerodynamic  Flows,  California  State 
University,  Long  Beach,  CA,  1981. 
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variable  with  discretisation  At  -  tn  -tn.  A  local  tine  linearisation  of 
requisite  formal  accuracy  is  introduced*  This  serves  to  define  a  linear 
differential  operator  L  such  that 


D{|"+1)  -  D(J")  ♦  Ln  (I"-1,1  -  ?*>  ♦  0(At2) 


(12) 


Thus,  Bq.  (11)  can  be  written  as  the  linear  system 

(I  -  0AtLn)  (J"+1  -  |")  -  AtD(|”)  ,  (13a) 

j"+1  -  .  I”*1  -  +  0(At3,  (P-1/2) At2)  ,  (13d) 


where  ^"+1  differs  from  Jn+1  by  the  linearisation  error. 

To  obtain  a  split  scheme,  the  multidimensional  operation  L  is  divided 
into  m  "one-dimensional”  suboperators  L«Lj+L2+.  .  .  «  1^,  which  are  usually 
associated  with  coordinate  directions.  For  simplicity,  only  two  suboperators 
are  considered  here i  Inclusion  of  additional  suboperators  is 
straightforward.  Bq.  (13)  can  then  be  replaced  by  the  approximate 
factorisation 


(I  -  PAtLj")  (I  PAtL  ")  <|*  -  I")  -  At  D($")  ,  (14a) 

-  I"  -  £n+1-  0[At3,  (P-1/2) At2)  ,  (14b) 


♦*  +n+1  f*  fn+1 

where  differs  from  £  by  the  factorization  error.  Although  $  ,  , 

and  I"  1  are  interchangeable  without  formal  loss  of  accuracy,  the  distinction 

is  worthwhile,  since  they  entail  different  types  of  error  which  in  practice 

nay  differ  widely  in  magnitude.  Note  that  the  approxiuate  factorization  (14a) 

does  not  represent  a  simplification  of  Bq.  (13)  until  it  is  written  in  a  split 

form. 


The  most  obvious  splitting  of  Bq.  (14a)  is 


(I  -  PAtl^”)  $  -  At  D(|")  ,  (15a) 

(I  -  0AtL2n)  (|*  -  J°)  -  $  ,  (15b) 

where  $  is  an  intermediate  quantity  whose ^physical  significance  follows  from 
its  definition  (15b),  which  implies  that  4*  approximates  the  time  increment 
|"+1  -  to  order  At.  If  spatial  derivatives  appearing  in  and  I<2  are 
replaced  by  three-point  difference  formulas,  than  each  step  in  Bqs.  (15a)  and 
(15b)  can  be  solved  by  a  block-tridiagonal  elimination. 
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The  derivation  o£  the  algorithm  is  incomplete,  however,  sin^e  Bq.  (15a) 
cannot  be  jolved  for  y  without  deriving  boundary  conditions  for  4  from  th^se 
given  for  4.  If  function  values  of  41  are  given^  then  boundary  values  of  4  can 
always  be  derived  from  the  definition  (15b)  of  41  for  use  in  Bq.  (15a).  In 
practice,  however,  more  complex  boundary  conditions,  such  as  normal 
derivatives,  may  be  specified;  and  so,  here  we  consider  a  much  more  general 
nonlinear  boundary  condition  which,  after  linearization  as  in  Bq.  (12),  can  be 
written  in  the  form  Bn($n+1)  -  gft,^1).  Here,  Bn  is  a  linearized  operator 
which  may  include  the  same  derivatives  as  L2n,  and  g  is  given.  Applying  the 
operator  B  to  Bq.  (15b)  gives 

n  +  n  ♦"i  n  n  -mi 

3  (4)  -  B  <♦  -  ♦  )  -  0At  b"[L  "(4  -  f1)]  .  (16) 

n  +*  +n  ~  n+1  n 

Since  B  (4  -  4  )  ■  9  -  g  without  formal  ^oss  of  accuracy,  Bq.  (16)  can  be 

usod  to  derive  exact  boundary  conditions  for  4  from  the  given  boundary 
conditions  provided  B°  and  1*^'  commute,  which  is  unfortunately  often  not  the 
case,  lhe  need  for  commutativity  occurs^because  L2n  4  cannot  be  computed  at 
this  step  in  the  algorithm,  whereas  Bn(4  )  can  be  replaced  by  the  given  gn+1 . 

When  Bn  and  L2n  do  not  commute,  an  approximate  solution  tjl  can  be 
computed  instead  of  $,  where  $  satisfies 

_n  +*  n+1  n  .  „  , 

B  (4  )  -  g  -  g  ,  (17a) 

-  4  +  0[At«4  -  |”l)  ,  (17b) 

instead  of  Bq.  (16).  Hiis  is,  of  course,  equivalent  to  using  uncorrected 
(i.e^,  "physical")  boundary  conditions  for  ($  -  $")  as  boundary  conditions 

for  4.  Substituting  Bq.  (17b)  into  Bqs.  (15)  shows  that  an  additional  error 
of  order  AtlJ  -  $"1  is  introduced  by  the  use  of  "uncorrected"  boundary 
conditions  as  in  Bq.  (17a).  Hie  coverall  accuracy  in  terms  of ^approximating 
Bq.  (11)  is  0[At^,  (0-1/2)  At,  1$  -  4°*^  where,  in  turn,  14  -  4°I  is  of 

order  At.  Note,  however,  that  in  the  steady  state  (4  -  4”)  2  0  '  and 

Bq.  (17a)  becomes  an  exact  boundary  condition.  Consequently,  leaving  boundary 
conditions  uncorrected  as  in  Bq.  (17a^fc  introduces  no  error  in  steady 
solutions,  which,  in  Jurn,  satisfy  D(4  )  *  4»  It  is  worth  noting  that  unless 
the  solution  4  of  D(4  )  <*  0  is  unique,  the  steady  solution  4  obtained  need 
not  be  the  same  as  would  be  obtained  by  repeating  Bq.  (11)  until  |n+1  -  411  ■ 
Hiis  completes  the  derivation  of  the  scheme  (15a)  and  (15b)  and  its  boundary 
conditions. 

Hiree-point  central  spatial  differences  have  been  employed  at  all 
interior  grid  points  in  the  difference  mesh.  At  the  boundaries  of  the 
computational  domain,  second-order  accurate  three-point,  one-sided,  spatial 
differences  are  used  to  represent  first  derivatives  where  required.  An 
artificial  dissipation  terra  based  upon  a  cell  Reynolds  number  criterion  has 
been  selectively  introduced  into  the  scheme  to  allow  calculations  to  be 
performed  at  high  local  values  of  the  cell  Reynolds  numbers. 


+0 


The  coordinate  system  for  interior  ballistios  calculations  must  have  the 
ability  to  enlarge  the  physical  extent  of  the  computational  domain  as  the 
projectile  moves  down  the  gun  tube.  To  accommodate  this  constraint ,  an 
accordion-type  mesh  is  used  in  the  axial  direction)  i.e.,  the  first  and  last 
axial  grid  points  are  attached  to  the  breech  and  projectile,  respectively)  and 
the  mesh  expands  as  the  projectile  accelerates  down  the  gun  tube.  This 
transformation  is  n  ■  (z-zQ)/(z  (t)-zQ),  where  *0,  z,  z  (t)  are  the  axial 
physical  distances  to  the  breech,  to  the  grid  point,  and  to  the  projectile 
base,  respectively.  Transformations  are  required  to  refine  the  computational 
mesh  to  regions  of  largo  gradients,  such  as  near  the  walls  and  where 
propagating  pressure  waves  appear  during  the  ignition  phase  of  the  ballistic 
cycle.  These  exponential  grid  point  concentrations  are  listed  in  Refs.  7  and 
8. 


IV.  THE  MODEL  PROBLEM 

Because  the  phenomena  occurring  in  a  two-phase,  multidimensional  interior 
ballistics  environment  are  extremely  complex,  an  idealized  gun  system,  a  one- 
phase  flow  in  the  Lagrange  gun,  has  been  simulated  to  facilitate  the 
understanding  of  some  of  the  basic  mechanisms  present  in  a  real  gun.  The 
Lagrange  gun  is  a  smooth  tube  of  constant  radius  which  is  closed  at  one  end  by 
the  breech.  The  combustion  chamber,  formed  by  the  breech,  the  initial 
position  of  the  flat-based  projectile,  and  the  tube  wall,  is  filled  with  a 
high  pressure,  high  temperature  gas  which  is  at  rest.  The  ballistic  cycle  of 
this  idealized  gun  resembles  that  of  a  real  gun  if  the  propellant  is  assumed 
to  be  completely  burned  before  the  projectile  moves.  The  projectile  is 
released  at  time  t-0 .  Frictional  forces  between  the  tube  wall  and  the 
projectile  are  neglected. 

At  the  present  stage  of  the  development,  we  report  on  the  results  of 
three  simulations i  the  laminar  flow-adiabatic  wall,  the  turbulent  flow- 
adiabatic  wall,  and  the  laminar  flow-isothermal  wall.  While  recognizing  that 
the  flow  in  the  Lagrange  gun  is  turbulent  and  heat  loss  through  the  walls  does 
affect  the  flow  field,  these  simulations  isolate  the  effects  of  those 
phenomena  and  can  help  in  the  validation  of  future  codes.  A  turbulent 
isothermal  computation  is  in  progress.  The  core  flow  under  laminar-adiabatic 
assumptions  has  been  verified  by  an  analytic  solution  and  the  other  cases 
display  the  expected  trends  in  the  flow  field.  These  benchmark  calculations 
assume  cylindrical  symmetry  and  solve  the  continuity  equation,  Eq.  ( 1 ) ,  the 
Navier-Stokes  equations,  Bqs.  (3),  and  the  energy  equation,  Eq.  (5),  with  the 
porosity  set  to  one.  The  single  species  gas  obeys  the  Noble-Abel  equation  of 
state,  Eq.  (38),  Sutherland's  laws  for  viscosity,  Eq.  (25),  and  thermal 
conductivity,  Eq.  (36),  and  Fourier's  law  of  heat  conduction,  Eq.  (34), 
with  ot*1.  For  the  turbulent  flow  calculations,  these  equations  are 
augmented13  by  the  turbulent  kinetic  energy  equation,  Eq.  (6),  and  the 
equations  for  turbulent  viscosity,  Eqs.  (27),  turbulent  thermal  conductivity. 


13B.E.  Launder  and  D.B.  Spalding ,  "The  Numrioal  Computation  of  Turbulent 
Flowe,"  Computer  Methods  in  Applied  Meohanioe  and  Engineering,  Vol.  3, 
pp.  289-189,  1974 . 
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Bq.  (37),  and  the  turbulent  heat  flux  vector,  Bq.  (35).  The  mixing 

length  1  was  set  proportional  to  the  distance  from  the  confining  surface  and  a 

van  Driest-type  of  damping  was  added  to  account  for  viscous  effects. 

The  following  boundary  conditions  are  imposed:  at  the  breech, 
projectile,  and  tube  wall,  a  no-slip  condition  is  maintained,  and  the  density 
is  determined  from  the  normal  momentum  equation.  The  boundary  conditions 
along  the  centerline  of  the  tube  are  the  symmetry  conditions 

u«0,  3w/3r-0,  3T/3r-0,  and  3p/3r>0.  The  temperature  boundary  condition  at  the 
solid  boundaries  is  that  the  normal  derivatives  of  the  temperature  are  set 
equal  to  zero  for  the  adiabatic  simulations  and  alternatively  the  temperature 
is  set  equal  to  300  K  for  the  isothermal  case.  In  the  turbulence  simulation, 
the  turbulent  kinetic  energy  and  the  mixing  length  are  zero  on  all  solid 
surfaces.  Tb  start  the  turbulence  calculations,  a  small  turbulent  kinetic 
energy  must  be  assumed  to  be  present  at  time  t-0.  The  gun  tube  geometry,  gas 
properties,  and  initial  conditions  are  the  same  for  each  simulation  and  ate 
given  in  Table  1 . 


TABLE  1.  PARAMETERS  FOR  THE  LAGRANGE  GUN  SIMULATIONS 


Bore  Diameter 

Combustion  Chamber  Length 

Maximum  Travel  of  Projectile 

Projectile  Mass 

Covolume 

Ratio  of  Specific  Heats 
Initial  Gas  Pressure 
Initial  Gas  Temperature 


20  mm 
0.175  m 
1.115  m 
0.120  kq 

1.08  x  10~3  m3A9 
1.271 
300  MPa 
3000  K 


The  axial  domain  is  uniformly  divided  with  49  mesh  lines.  The  radial 
domain  is  nonun if ormly  divided  with  19  mesh  lines  and  with  a  mesh  line 
concentration  near  the  tube  wall,  where  the  first  mesh  line  is  at  the  wall, 
the  second  mesh  line  is  7.7  pm  away  from  the  wall,  the  tenth  mesh  line  is  0.8 
mm  away  from  the  wall,  and  the  nineteenth  mesh  line  is  the  centerline  of  the 
gun  tube. 


A.  Laminar  Flow-Adiabatic  Wall  Simulation 


The  dominant  features  of  the  one-phase  flow  within  the  Lagrange  gun  are 
the  existence  of  the  rarefaction  wave  propagation  and  the  development  of  the 
velocity  and  thermal  boundary  layers  along  the  gun  tube  wall.  The  rarefaction 


16 


'j#av«  is  generated  by  tha  notion  of  tha  projaotila  and  nay  travaraa  tha 
distance  between  tha  breech  and  projectile  one  or  more  tines  before  the 
projectile  exits  the  tuba.  The  ALPHA  confuted  pressure  histories  at  the 
center  of  the  breech  and  projectilo  base  are  given  in  Figure  2.  As  predicted 
in  1793  by  Lagrange'4,  using  many  simplifying  assumptions  awd  a  lumped 
parameter  analysis,  the  pressure  at  the  breech  is  always  larger  than  that  at 
the  projectile  base  for  this  idealised  gun  system.  The  position  of  the 
rarefaction  wave  is  indicated  by  a  very  rapid  change  in  the  slope  of  the 
pressure  curve.  In  numerical  solutions,  such  changes  are  often  smeared  out. 
The  first  arrival  of  the  rarefaction  wave  at  the  breech  and  projectile  base  is 
captured  by  the  numerics  as  sasn  in  Figure  2.  In  the  core  region,  the  flow  is 
predominantly  axial  and  at  early  tines  is  nearly  isentropic.  A  ona- 
dimensional  analytic  solution  for  the  isentropic  case  ia  developed  in 
Ref.  15.  The  computed  pressure  values  at  the  time  of  the  rarefaction  wave  at 
the  breech  end  projectile  are  compared  to  the  analytically  determined  values 
in  Figure  2.  The  differences  between  the  values  are  2.8  percent  at  the  breech 
and  3.4  percent  at  the  projectile  base.  The  numerical  results  were  obtained 
with  a  constant  time-step  of  10  vs. 


TIME  (mi) 

Figure  2.  Pressure  Histories  at  the  Breech( - )  and  at  Projectile 

B*$e( — ),  and  Their  Comparisons  with  Analytically  Predicted 
Initial  Arrival  Times  of  the  Rarefaction  Wave 


Comer ,  Theory  of  Interior  Ballietioe  of  Gune,  Wiley ,  New  York, 
pp.  339-256,  1950. 

1SE.H.  Love  and  F.  B.  Pidduok,  " Lagrange's  Ballistic  Problem,  "  Phil.  Trane. 
ROU .  Soc.,  Vol.  222,  pp.  167-226,  1921-22 . 
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The  gun  gnome try,  gas  parameters,  initial  conditions,  and  boundary 
conditions  used  in  the  simulation  are  identical  to  those  used  by  Heiser  and 


Hensel 


1i 


.1  *! 


Despite  these  similarities,  the  results  of  the  simulations  differed 


greatly’'.  Subsequent  results  from  West  Germany^  show  agreement  between  the 
simulations.  For  example,  at  a  muzzle-time  of  2.6  ms,  the  differences  between 
these  calculations  at  the  projectile  base  are  0.4  percent  in  the  temperature, 
0.8  percent  in  the  axial  displacement  values,  and  2.2  percent  in  the  velocity 
values . 


The  other  dominant  features  of  the  flow  field  within  the  Lagrange  gun  are 
the  boundary  layers  along  the  gun  tube  wall .  The  velocity  boundary  layer 
profiles  at  time  2.4  ms  are  given  in  Figure  3,  where  z and  w  denote  the 
axial  position  and  velocity  of  the  projectile,  respectively.  The  99  percent 
velocity  boundary  layer  thickness  has  a  maximum  value  along  the  tube  of  0.19 
mn.  A  more  physically  meaningful  measurement  for  the  boundary  layer  thickness 
is  the  displacement  thickness  6.  For  compressible  flows  with  cylindrical 
geometry,  it  is  defined  at  time  t  and  position  z  as 


R 

6(z,t)  *  /  [i 
o 


w(r,z,t)  p(r,z,t) I  £ 

w  (z,t)  p  (z,t)  J  R 
c  c 


(18) 


where  the  subscript  c  denotes  the  value  at  the  centerline,  w  is  the  axial 
velocity,  p  is  the  density,  and  R  is  the  tube  radius.  The  displacement 
thickness  values  given  in  Ihble  2  are  approximately  27  percent  of  the  99 
percent  velocity  boundary  layer  thickness  values.  At  2.4  ms,  the  maximum 
displacement  thickness  occupies  less  than  0.6  percent  of  the  gun  diameter. 

The  thermal  boundary  layers  at  2.4  ms  are  given  in  Figure  4.  The  core 
temperatures  decrease  as  a  function  of  the  distance  from  the  breech  due  to  the 
motion  of  the  projectile.  The  temperature  profile  at  the  projectile  base  is 
approximately  constant  in  the  radial  direction.  Along  the  adiabatic  tube 
wall,  the  viscous  heating  raises  the  wall  temperature  to  a  maximum  value  away 
from  the  projectile  base,  e.g.,  in  Figure  4  at  Z  ■  0.805  m.  The  wall 
temperature  then  decreases  from  this  concentrated  region  of  high  temperature 
as  one  proceeds  to  the  breech. 


7  ft 

R.  Heiser  and  D.  Heneel ,  " Calculation  of  ti le  Axi symmetric  Unsteady 
Compressible  Boundary  Layer  Flow  Behind  a  Moving  Projectile, "  Proceedings  of 
tiie-Fourth  International  Symposium  on  Ballistics.  1978,  American  Defense 
Preparedness  Association,  Washington,  DC . 

J.A.  Schmitt  and  T.L .  Mann ,  " Calculation  of  the  Compressible  Flow  in  the 
Lagrange  Gun  by  the  Interior  Ballistics  Algorithm,  ALPHA."  Proceedings  of 
the  DEA-G-10  60  Meeting  at  Ealin  Air  Force  Base.  October  1980,  Air  Force 
Armament  Division,  Eglin  Air  Force  Base ,  Florida,  pp.  4-25. 

7  ft 

R.  Heiser  and  D.  Hensel,  " Beredhnung  der  Gasstromung  in  einem  Waffenrdhr  mit 
Hilfe  dee  Zweidimensionalen  AMI-Modelle ,"  El/81,  January  1981 ,  Emst-Mach- 
Inetitut,  Abteilung  fuer  Ballistik,  Weil  am  ftiein.  West  Germany. 
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DISTANCE  FROM  BREECH  (m) 


Figure  3.  Velocity  Boundary  Layer  Profiles  for  the  Laminar 
Flow-Adiabatic  Wall  Calculation 


TABLE  2.  DISPLACEMENT  THICKNESSES  IN  THE  LAGRANGE  GUN 

AT  TIME  t«2.4  as 

Displacement  Thickness  (mm) 


Distance 
from  Breech 
(m) 

Laminar- 

Adiabatic 

Case 

Turbulent- 

Adiabatic 

Case 

Laminar- 

Isothermal 

Case 

.073 

.052 

.094 

.005 

.268 

.052 

.161 

.002 

.439 

.052 

.193 

.002 

634 


050 


210 


002 


Figure  4.  Thermal  Boundary  Layer  Profiles  for  the  Laminar 
Flow-Adiabatic  Wall  Calculation 


B.  Turbulent  Flow- Adiabatic  Wall  Simulation 


Indications  are  that  the  flow  In  a  gun  barrel,  even  In  an  Idealized 
Lagrange  gun,  Is  turbulent.  Thus,  a  k-Jt  turbulence  model  was  Incorporated 
Into  the  computational  scheme  and  the  results  compared  with  the  laminar  runs 
to  determine  If  there  are  significant  differences. 

Perhaps  the  most  interesting  result  of  the  turbulence  simulation  is  the 
displacement  layer  (Table  2).  At  muzzle  time,  at  a  distance  of  0.8  m  from  the 
breech,  the  displacement  thickness  is  0.2  mm.  The  velocity  profiles. 

Figure  5,  are  typical  of  turbulent  flow.  As  a  conse*'  nee  of  the  viscous 
heating  along  the  adiabatic  wall,  there  is  a  negatlv  emperature  gradient  in 

the  direction  away  from  the  wall;  for  example,  at  0.8  m  from  the  breech,  the 
wall  temperature  is  1742  K  while  the  centerline  is  still  at  1690  K  (Figure 
6).  Note  that  the  scale  of  the  ordinate  in  the  figure  is  so  large  that  the 
actual  grid  spacing  of  the  computational  domain  can  no  longer  be  discerned 
near  the  wall.  An  enlarged  view  would  show  that  the  adiabatic  wall  condition 
3T 

■  0  is  strictly  maintained.  The  level  of  turbulence,  measured  by  the 

turbulent  kinetic  energy,  increases  from  the  breech  toward  the  projectile. 

The  mixing  action  of  the  turbulence  transports  heat  from  the  viscous  region 
near  the  wall  into  the  core  flow  which  is  reflected  in  the  temperature 
gradient  seen  in  Figure  6. 
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Figure  5*  Velocity  Boundary  Layer  Profiles  for  the  Turbulent 
Flow-Adiabatic  Wall  Calculation 


Figure  6.  Thermal  Boundary  Layer  Profiles  for  the  Turbulent 
Flow- Adiabatic  Wall  Calculation 
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Laminar  Flow- Isothermal  (Cold  Wall)  Simulation 


The  simulation  of  the  laminar  flow  in  the  Lagrange  gun,  under  isothermal 
wall  conditions,  determines  the  effects  of  heat  losses  on  the  dynamics  of  the 
flow  field  and  the  motion  of  the  projectile.  This  simulation  is  a  reasonable 
approximation  of  the  heat  transfer  for  the  real  case  of  a  turbulent  flow  with 
heat  conducting  walls.  The  calculations  predict  that  approximately  20  percent 
of  the  total  energy  available  in  the  gas  is  lost  through  the  cold  walls. 

These  results  are  in  agreement  with  predicted  heat  losses  from  conventional 
propellants  in  medium  caliber  guns1. 

Radial  dependence  of  the  axial  velocity  profiles  as  a  function  of 
distance  from  the  breech  of  the  gun  is  shown  in  Figure  7.  The  maximum  value 
of  the  99  percent  velocity  boundary  layer  thickness  is  calculated  to  be  0.15 
mm.  The  corresponding  displacement  thicknesses  are  given  in  Table  2.  The  99 
percent  velocity  boundary  layer  thicknesses  of  this  case  are  comparable  to 
those  of  the  laminar  flow-adiabatic  wall  simulation.  However,  the 
corresponding  displacement  thicknesses  differ  by  over  an  order  of  magnitude 
because  the  density  variation  is  taken  into  account  in  the  displacement 
thickness  calculation  whereas  it  is  not  in  the  99  percent  velocity  boundary 
layer  computation.  Se^  Eq.  (18).  The  gas  density  near  the  wall  for  the 
isothermal  cold  wall  case  is  significantly  different  from  the  adiabatic  case 
because  the  pressure  values  are  approximately  the  same  but  the  temperature 
values  differ  greatly,  in  general,  the  thermal  boundary  layer  in  the  cold 
wall  simulation  remain  extremely  steep  over  the  time  span  of  in-bore  travel 
(Figure  8) . 


Figure  7.  Velocity  Boundary  Layer  Profiles  for  the  Laminar 
Flow- Isothermal  (Cold  Wall)  Calculations 


Figure  8.  Thermal  Boundary  Layer  Profiles  for  the  Laminar 
Flow- Isothermal  (Cold  Mall)  Calculation 


During  the  early  transients,  the  discontinuity  In  the  Initial  conditions 
for  the  temperature  creates  a  strong  transient  "thermal  discontinuity"  In  the 
chamber  cavity.  This  results  from  the  drop  In  temperature  through  the  thin 
layer  of  gas  adjacent  to  the  wall,  which  Is  needed  to  satisfy  the  boundary 
condition  while  keeping  the  pressure  at  300  MPa.  Consequently,  severe 
gradients  in  the  gas  density  near  the  walls  occur  and  the  convection  and 
diffusion  of  the  mass  to  these  regions  persist  for  most  of  the  ln-bore  travel 
time. 


D.  Comparison  Among  the  Simulations 

From  the  Results  we  see  that  there  are  some  significant  differences  In 
the  flow  fields  of  the  three  cases.  Calculations  with  a  turbulence  model 
reveal  marked  departures  from  the  laminar  case  in  the  details  of  the  flow, 
although  the  overall  trends  are  similar.  The  isothermal  case  differs  on  both 
the  micro  and  the  macro  scale.  For  the  adiabatic  wall,  the  velocity  profile 
normal  to  the  tube  wall  is  more  gradual  for  a  turbulent  calculation  than  the 
corresponding  laminar  one  with  a  corresponding  thickening  of  the  displacement 
thickness.  In  fact,  for  example,  at  a  distance  of  0.4  m  from  the  breech 
(Figure  9)  the  turbulent  displacement  thickness  is  0.19  mm  while  the 
corresponding  laminar  value  is  0.05  mm.  Due  to  the  eddy  motion  in  the 
turbulence,  there  is  a  considerable  mixing  of  fluid  leading  to  a  lower 
temperature  at  the  wall  and  a  more  rapid  transport  of  heat  into  the  core 
(Figure  4  and  6).  In  the  cold  wall  case,  the  displacement  thlcknes3  Is  very 
thin,  of  the  order  of  0.005  om  at  10  on  from  the  breech  near  projectile  exit 
time,  and  gets  progressively  thinner  as  one  moves  toward  the  projectile 
(Figure  9). 
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Figure  9.  Comparisons  of  the  Displacement  Thicknesses  for  the  Three 
Simulations  of  the  Model  Lagrange  Gun  Problem  at  t=2.4  ms 


From  Figure  10,  we  see  that  the  projectile  velocity  in  the  turbulent  case 
is  slightly  lower  than  in  the  laminar  case  and  the  velocity  difference 
increases  as  the  projectile  travels  down  the  tube,  confirming  the  fact  that 
energy  is  taken  up  by  the  turbulent  motion  and,  thus,  demonstrating  the 
consistency  of  the  results,  the  difference  between  the  adiabatic  and  the 
isothermal  cold  wall  projectile  velocities  for  the  laminar  simulation  is  quite 
large  and  increases  with  time.  Due  to  the  loss  of  energy  through  the  solid 
surfaces  during  in-bore  travel  time,  the  muzzle  velocity  in  the  cold  wall  case 
is  50  m/s  less  than  in  the  adiabatic  case.  Furthermore,  this  heat  loss  causes 
a  notable  temperature  difference  along  the  centerline  between  these  cases 
which  can  be  seen  by  comparing  Figures  4  and  8.  At  about  70  percent  of  the 
distance  from  the  breech  to  the  projectile  base,  the  temperature  difference  is 
approximately  110  K.  In  the  adiabatic  cases  (Figures  4  and  6)*,  heat  generated 
by  the  viscous  heating  at  the  wall  is  transported  and  conducted  into  the  flow 
leading  to  the  positive  temperature  gradient  as  one  approaches  the  wall. 
Whereas  in  the  isothermal  case  (Figure  8),  heat  is  lost  through  the  wall 
leading  to  a  negative  temperature  gradient. 


V.  CONCLUSIONS  AND  FUTURE  EFFORTS 

An  advancement  has  been  made  in  the  state-of-the-art  in  modeling 
computationally  the  multidimensional  Interior  ballistics  of  high-velocity 
guns.  Hie  numerical  LBI  scheme  is  used  to  solve  the  axisymmet.ric,  two-phase, 
averaged  Navier-Stokes  equations  which  describe  the  transient  flow  behind  an 
accelerating  projectile.  However,  in  the  first  application  of  this  algorithm, 
the  interior  ballistics  flow  has  been  idealized  to  a  single-phase  case;  and 
three  simulations  are  presented:  a  laminar  flow  with  an  adiabatic  wall 


TIME  (ms) 

Figure  10.  Differences  In  Projectile  Velocities  as  a  Function  of 
Time,  Flow,  and  Boundary  Conditions 


condition,  a  turbulent  flow  with  an  adiabatic  wall  condition,  and  a  laminar 
flow  with  an  isothermal  cold  wall  condition.  Despite  a  relatively  coarse 
mesh,  the  rarefaction  wave  has  been  captured;  the  projectile  motion  has  been 
successfully  coupled  to  the  two-dimensional  flow  field;  and  the  velocity  and 
thermal  profiles,  as  well  as  displacement  thicknesses,  have  been  determined 
for  all  the  cases.  Comparisons  with  analytical  results  and  others,  when 
available,  ha  ye  resulted  in  excellent  agreement.  Furthermore,  confidence  in 
these  results  were  enhanced  through  time-step  and  mesh  refinements. 

Past  attempts  to  model  the  viscous  flow  within  an  idealized  ballistic 
cycle  has  been  confined  to  boundary  layer  studies  via  the  boundary  layer 
equations1’.  They  predict  trends  analogous  to  those  reported  here;  such  as, 
the  thickening  of  the  boundary  layer  along  the  tube  wall  under  turbulent 
conditions.  However,  because  those  approaches  fail  to  consider  the  entire 
flow  and  the  flow-projectile  interaction,  they  cannot  simulate  the  decrease  in 
the  projectile  velocity  due  to  turbulence  or  heat  transfer  to  the  wall.  These 
effects  can  be  significant. 

In  the  ono-phase  simulations  reported  here,  important  aspects  of  the  two- 
phase  flow  were  excluded  from  the  initial  validation  of  the  basic  models  and 
the  algorithm.  The  next  step  in  the  development  of  the  viscous  simulation  of 
the  ballistic  cycle  is  the  inclusion  of  the  solid  phase  in  the  Lagrange  gun 
and  the  validation  of  the  corresponding  two-phase  models.  To  this  end, 
benchmark  experiments  in  a  well-controlled  ballistic  environment  are  being 
designed  in  conjunction  with  the  modeling  efforts.  TVo  of  the  most  pressing 
needs  in  the  modeling  effort  in  two-phase  flow  are  the  development  of  models 
of  turbulence  for  the  second  phase  and  a  more  realistic  description  of  the 
interaction  between  phases . 


^E.P.  Bartlett,  L.W.  Andereon ,  and  R.M .  Kendall,  "Time-Dependent  Boundary 
Layere  uri-th  Application  to  Gun  Barrel  Heat  Transfer."  Proceedings  of  Tuelfth 
Heat  Transfer  and  Fluid  Mechanics  Institute .  Stanford  University  Preea, 
pp.  262-278,  1972. 
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In  this  Appendix,  we  list  the  constitutive  laws  and  the  correlations 
needed  to  close  the  two-phase  model.  The  mass  source  due  to  the  propellant 
burning  is 
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where  the  grain  surface  area,  grain  volume,  and  regression  rate  are 
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respectively.  The  constants  ,  Bj,  and  n  are  known  for  a  given  propellant. 
The  gas  phase  stress  tensor,  assuming  a  Newtonian  fluid,  is 


w  «  2pD  -  (-|p  -  1^)  V  •  U  1!  , 


(23) 


where  p  is  the  molecular  viscosity,  Kfi  is  the  bulk  viscosity  and  II  is  the 
identity  tensor.  The  rate  of  strain  tensor  B  is  given  as 
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The  molecular  viscosity  p  is  determined  from  Sutherland's  law 
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where  S1  *  110  K  and  the  subscript  zero  indicates  some  reference  quantity. 
The  bulk  viscosity  Kg  is  assumed  to  be  zero.  The  turbulent  stress  tensor  is 
modeled  using  an  isotropic  eddy  viscosity  formulation 
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where  j l  is  the  mixing  length  and  c  ie  a  function  of  the  turbulent  Reynolds 
number.  The  interphase  drag  relation  <F>  is  given  by 
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where  a  is  the  settling  porosity  of  the  propellant  bed.  The  intergranular 
stress  relation  R  which  is  independent  of  the  loading  density  history  is 
gi v  an  by 
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else 


where  a  is  the  speed  of  sound  in  the  solid  phase  which  is  specified  for  a 
given  propellant.  The  dissipation  function  ♦  is  given  by 
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The  interfacial  energy  transfer  A  is  given  by 
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where  hCQ  is  the  energy  released  per  unit  mass  due  to  combustion  of  the 
propellant.  The  interfacial  heat  transfer  <q>  is  given  by 
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where  e  «  propellant  emissivity,  a  *  Stefan-Boltzmann  constant,  and 
Pr  *  yc^/K.  The  surface  temperature  of  the  propellant  TpS  is  computed  from 
the  solution  of  0q.  (9).  The  laminar  and  turbulent  heat  flux  vector  can  be 
modeled  as 


*►  voi 

q  -  -  K  [VT  - -  (T^-T) ] 


(34) 


and 


+T  T  ,  _  Va  .  , 

q  “  -  K  [VT - -  ( T^-T)  ] 


(3S) 


respectively,  where  »  0.5(T+T  ).  The  thermal  conductivity  K  is  determined 

from  Sutherland's  law 
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where  S2  ■  194  K  and  subscript  zero  indicates  some  reference  quantity.  The 
turbulent  thermal  conductivity  is  given  by 
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where  Prejf  *  0.9.  The  Noble-Abel  equation  of  state  is 
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where  h  is  the  covolume  factor  which  provides  a  correction  to  the  perfect  gas 
equation  of  state  needed  for  gases  with  large  density  and  is  the  universal 
gas  constant.  Other  thermodynamic  properties  are 
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=  speed  of  sound  in  the  propellant,  m/s 

-  linearized  operator 

»  specific  heat  at  constant  pressure,  J/(kg»K) 

»*  specific  heat  at  constant  volume,  J/(kg»K) 

«  propellant  thermal  diffusivity,  m  /s 
=  regression  rate  of  solid  phase,  m/s 
=  spatial  differential  operator 

■  rate  of  strain  tensor,  1/s 

=  interphase  drag  per  unit  area  of  solid  phase.  Pa 

■  specific  enthalpy,  J/kg 

-  identity  matrix 

■  specific  turbulent  kinetic  energy,  J/kg 
«  bulk  viscosity,  Pa»s 

=»  laminar  thermal  conductivity,  W/(m*K) 

*  turbulent  mixing  length,  m 

=  linearized  differential  operator 

a>  reciprocal  molecular  weight  of  gas  mixture,  mol/k9 
=  pressure.  Pa 

=  laminar  heat  flux  vector,  J/(m  »s) 

■  turbulent  heat  flux  vector,  J/(m  *s) 

*  interfacial  heat  transfer,  W/(m  ) 

«  isotropic  intergranular  stress.  Pa 
*»  universal  gas  constant,  J/(mol*K) 

=  radial  distance  from  the  centerline  in  the  gun  tube,  m 
=  propellant  radius,  m 

2 

»  propellant  surface  area,  m 

-  0.9 

■  time,  s 

*  gas  temperature,  K 
=  propellant  temperature,  K 
=  radial  component  of  gas  velocity,  m/s 
=  gas  velocity  vector,  m/s 
=  propellant  velocity  vector,  m/s 
=  propellant  volume,  m^ 

=*  axial  component  of  gas  velocity,  m/s 
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